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ABSTRACT 

We show that collisions with stellar-mass black holes can part ially explain t he ab - 
sence of bright giant stars in the Galactic Centre, first noted bv lGenzel et alj (|1996f) . 
We show that the missing objects are low-mass giants and AGB stars in the range 
1 — 3M©. Using detailed stellar evolution calculations, we find that to prevent these 
objects from evolving to become visible in the depleted K bands, we require that they 
suffer collisions on the red giant branch, and we calculate the fractional envelope mass 
losses required. Using a combination of Smoothed Particle Hydrodynamic calculations, 
restricted three-body analysis and Monte Carlo simulations, we compute the expected 
collision rates between giants and black holes, and between giants and main-sequence 
stars in the Galactic Centre. We show that collisions can plausibly explain the missing 
giants in the 10.5 < K < 12 band. However, depleting the brighter (K < 10.5) objects 
out to the required radius would require a large population of black hole impactors 
which would in turn deplete the 10.5 < K < 12 giants in a region much larger than 
is observed. We conclude that collisions with stellar-mass black holes cannot account 
for the depletion of the very brightest giants, and we use our results to place limits on 
the population of stellar-mass black holes in the Galactic Centre. 
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> . 1 INTRODUCTION 



In common with many other galaxies, the Milky Way hosts 
at its centre a supermassive black hole, Sagittarius A* , 
whose mass is » 4 x 10 6 M Q (e.g. ISchodel et ail 120031 ). 
Surrounding Sgr A* is the Galactic Centre cluster, the 
study of which offers unique insights into many areas of 
astrophysics, such as the dynamics of star clusters with 
central massive objects, the formation history of supermas- 
sive black holes and the emission of gravitational waves by 
Extreme Mass-Ratio Inspirals (for a recent review of this 
topic, see lAmaro-Seoane et alj (|2007h . 

The study of the Galactic Centre cluster is hindered 
by ~ 30 magnitudes of intervening visual extinction. Nev- 
ertheless, a great deal has been learned about the stellar 
population in the central few parsecs by observations made 
in X-rays and in the near infra-red, in which the extinction 
is only ~ 3 magnitudes. 

The Galactic Cent re Cluster is one of the densest 
stellar systems known. (|Genzel et alj fl996l ) concluded that 
the Cluster had a dense core with a radius of w 0.38pc 
and a densit y of 4 x 10 6 Mp)pc~ 3 . Late r meas urements by 
iGenzel et all (|2003l ) and ISchodel et all (|2007l ) found that 
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the cluster was more cusp-like and reported still higher 
number densities, reaching « 10 7 pc -3 at a radius of O.lpc. 
These number densities are well in excess of those of globu- 
lar clusters (< lO 6 M pc" 3 , lHarrislll996h or young cluster s 



such as the Arches (3 x 10 5 MmPC~ 3 . iFiger etahl Il999h 
The c usp is composed lar g ely of old and interm ediate-age 
stars l|Genzel et all 120031 ; ISchodel et ail 120071 ). However, 
mixed with this population a re two unusual groups of 
young stars. The S-stars fe.g. iGhez et al.ll2005l ). orbiting 
within ~ 0.01 pc have attracted considerable attention 
because they appear to be young B-stars living in a volume 
in which star formation should be strongly suppressed by 
the tidal field of Sgr A*. Alternative suggestions for the 
origin of the S-stars include the inward migration and 
dissolution of a dense stellar cluster fe.g. iGerhard 12001 



Hansen fc Milosavlievid 120031 ; iPortegies Zwart et al 



2003 



Giirkan fc Rasioll2005l) or the disruption of a mass ive binary 



injected from larger radii (|Gould fc Quiller]|2003l ). in which 
case they are young b ut formed elsewhere, and the tidal 
stripping of giant stars l|Davies fc K ing 2005), in which case 
the S-stars are locally formed but not young. 

Further out, between radii of ~ 0.04 and ~ 0.3 pc, 
there appear to be two discs of young stars oriented at 
large angles both to each other and to t he Galactic plane 
l|Genzel et al.ll2003l ; iPaumard et""ai1l2006l ). There is strong 
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evidence that these discs contain large nu mbers of massive 
stars and have an unusually flat IMF (|Paumard et al.l 
l200d ). 

In this paper, we investigate a peculiar property of 
the background cluster of older stars; within ~ 0.2pc, the 
Galactic Centre cluster is defici ent in bright giant stars. 
This was first documented by iGenzel et al.l |l996) who 
observed the Galactic Centre in the near infrared K-band. 
They used the presence or absence of CO bandhead 
absorption to distinguish late-type from early-type stars 
and divided their sample into three K-bands; 12 < K < 15, 
10.5 < K < 12 and K < 10.5. While the surface density 
of late-type stars in the faintest of these bands appears to 
continue smoothly inwards to very small radii, there are 
clear holes in the distributions of late-type stars in the 
other bands, with projected radii of ~ 0.08 and ~ 0.2 pc 
respectively. 

Since the Galactic C entre is a densely-populated re- 
gion, IGenzel et al.l 0996) proposed that the missing giants 
had been destroyed by collisions wit h main-sequence stars. 
This scenario was inv estigated by lAlexand"erl (|l999l ) and 
iBailev fc Davies] (1 19991) but th ese authors reached different 
conclusions. Alexander! l|l999l ) found that collisions with 



MS s tars could expla i n the observed depletion of giants, 
while IBailev fc Daviesl (l999) concluded that collisions with 
main sequence stars, white dwarfs and neutron stars would 
not do sufficient damage to a giant to prevent it becoming 
visible in the depleted K-bands. Their different conclusions 
largely arise from the criteria they used to determine 
whether or not an impact had done sufficient damage to 
prevent a given giant ap pearing in the de p leted K-bands 
in the Galactic Centre. IBailev fc Daviesl l|l999l ) studied 
encounters involving two very late giants. They found 
that the mass losses achieved by collisions with MS and 
compact-object impactors were small, no greater than 10% 
of the envelope. They also considered that removal of 'most 
if not nearly all of the envelope' was required to prevent 
the giant becoming visible in the depleted K-bands (as we 
will show later, this is largely correct for giants towards 
the tip of the giant branch). They therefore concluded 
that collisions with MS stars we re unable t o expl ain the 
depletion of giants. Conversely, lAlexanderl (|l999T l made 
the assumption that a collision with a single MS impactor 
would result in the total destruction of a giant's envelope if 
the ratio of the impact parameter to the giant radius, x e , 
were less than 0.25, and that an impact with a binary would 
destroy the giant envelope if x e < 1. Essentially takin g 
collisions t o be m ore effective than IBailev fc Daviesl l|l999T ) , 
lAlexanderl l| 19991 ) found that encounters could explain 
the depletion of giants within a radius of ~ 2 arcsec. By 
performing hydrodynamic modelling of encounters between 
MS impact ors and younger, l ess ex tended giants than those 
studied by IBailev fc Daviesl l|l999l) . we improve on and ex - 
tcnd the work of both of these authors. iDavies et al.l (|l998l ) 
investigated the possibility that collisions with binary MS 
stars or neutron stars might be responsible for removing 
the giants. They concluded that such encounters would 
not have a significant effect on the giant population in the 
Galactic Centre unless an unrealistically large proportion 
of binaries was assumed. 

On theoretical grounds, a large population of stellar- 
mass black holes is expected to exist in the Galactic 



Centre (e.g iMorris! 1 19931 : iMiralda-Escude fc Gould! l200d : 
iFreitag et ail 12006) (for re views of stella r rela xation 
processes around SMBH, see lAlexanderl r2006l . 120071 ). This 
population is of great interest, but is by its nature extremely 
difficult to study. In this paper, we explore the possibility 
of using the effects of collisions involving these black holes 
to study this population indire ctly. 

We extend the work of IBailev fc Daviesl (|l999l ) by 
considering collisions between giants and stellar-mass black 
holes, and also by allowing for multiple collisions of a single 
giant with black holes, MS stars or a combination of the 
two types of object. We then consider what can be learned 
from our results about the population of stellar mass black 
holes that surely exists in the Galactic Centre. 

In Section 2, we discuss the properties of the Galactic 
Centre stellar population and identify the stars that are 
missing. In Section 3, we derive the expected collision rates 
of giant stars with main-sequence stars and black holes. 
We describe our numerical methods in Section 4 and in 
Section 5 we explain how we have modeled the effects of 
mass loss on giant stars. In Section 6, we present the results 
of our simulations of collisions of giants with black holes 
and main-sequence stars. We discuss our results in Section 
7 and draw our conclusions in Section 8. 



MISSING STARS IN THE GALACTIC 
CENTRE 



Observations by IGenzel et alj (|l996T ) in three different 
infra-red K-bands showed that there is a region very close 
to Sgr A* devoid of bright giant stars. They detected no 
late-type stars with K magnitudes brighter than K — 10.5 
within w 0.2 pc and no late-type stars with K magnitudes 
in the range 12 > K > 10.5 within « 0.08 pc of the 
supermassive black hole. In a fainter band, 15 > K > 12, 
they observed that the surface density of late-type stars 
continued smoothly in to much smaller radii. For the 
remainder of the paper , we w ill refer to the three bands 
defined in IGenzel et al.1 (|l996Tl as the bright, midd l e and 
faint bands respectively. Foll owing Bailey fc Daviesl (|l999l ) 
and referring to Figure 6 in IGenzel et al.1 i| 19961 ), we esti- 
mate that the number of objects missing from the bright 
band is between eleven (assuming that the surface density 
distribution of these objects is flat inside 6 arcsec) and 
twenty-five (extrapolating the surface density inwards). 
Similarly, we estimate that the number of objects missing 
from the middle band is between six and eleven. The 
sharp drops in surface density at ~ 6 arcsec in the bright 
band and ~ 3 arcsec in the middle band argue against the 
depletion being caused by peculiarities in the IMF and 
we do not think it can be due to tidal stripping, as will 
be explained in a later paper. Instead, the nature of the 
depletion strongly implies a collisional origin. 



2.1 Identifying the missing stars 

In Figure [T] we show an HR diagram depicting the 
evolutionary tracks of 1,2,3,4, 8,15,20 and 30Mp) stars 
derived using the STARS code l|Eggletonl fl97ll : iPols et all 
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Figure 1. Hcrtzprung-Russell diagram showing tracks of 
1,2,3,4, 8, 15,20, 30 M Q stars with contours of K=10.5, 12, 15 
(corresponding to the boundaries of the bands observed in 
Gcn zel et afl dl996)l using main-sequence colours and bolometric 
corrections (dotted lines) and giant colou rs and bolometric correc- 
tions (dashe4_lines)from [Johnson (1966). A distance modulus of 
14.6 llGenzel et~al Jl996h and an extinction Ak of 3.0 l lRieke et al.l 
1 1989b have been assumed. Tracks are for solar-metallicity stars. 
Objects with K < 10.5 are depleted within ~ 0.2 pc, objects with 
10.5 < K < 12 are depleted within ~ 0.08 pc and objects with 
12 < K < 15 are not depleted. 



1 19951 ). We overlay K contours of 10.5, 12 and 15 magni- 
tudes calculated __using colours and bolometric corrections 
from [johnson ( 196fJ) and adopting a distance modulus 
of 14.6 llGenzel et al.l 1 19961) and an extinction Ak of 3.0 
l|Rieke et al.lll989l ) and assuming solar metallicity. Dotted 
lines show contours appropriate for main-sequence stars 
and dashed lines are contours appropriate for giants. 
This plot shows that objects brighter than K — 10.5 are 
low-mass AGB stars or high-mass red supergiants, while 
those whose K magnitudes lie between 12 and 10.5 are likely 
to be low-mass red giants and AGB stars or high-mass 
main-sequence (MS) stars. 

To obtain a more detailed picture of which segment 
of the Galactic Centre stellar population is missing, we 
calculate the length of time spent by each star within a 
given K-band. If we assume that the stellar population in 
the Galactic Centre has been forming stars continuously 
at a steady rate following a Miller-Scalo IMF for 14 Gyr, 
we can estimate the relative proportions of objects in 
a population of this age which are visible at a given K 
magnitude. We also attempted realisations using different 
assumptions. We found that a Salpeter IMF was unable 
to reproduce the numbers of giants that are observed 
outside the deplet ed regions in the Galactic Centre (e.g. 
iGenzel et all [l996) in the K-bands using a realistic total 
mass for the Galactic Centre cluster - the result is a giant 
population in which the bright band is overpopulated with 
respect to the middle and faint bands. We were unable to 
significantly improve the fit to the observed numbers of 



giants by varying the slope of the Salpeter mass function 
from its canonical value of 2.35 and the overabundance 
of bright giants becomes significantly worse for slopes 
flatter than 1.5. However, we note that, in accord with the 
Miller-Scalo IMF, Salpeter-like IMFs result in the bright 
K-band being largely populated by 2-4M0 objects and 
the middle band by 1-2M© objects, unless the slope was 
< 1.5. In addition, we found that assuming a burst of star 
formation (i.e. that the Galactic Centre stars are all the 
same age) was unable to reproduce the observed numbers 
of giants regardless of the time at which the burst occurred, 
since giants of a very narrow age (and therefore mass) range 
all have very similar K-band magnitudes. In particular, a 
single burst occurring > 10 Gyr ago would result in only 
giants of « IMq being visible now, leaving the bright band 
almost unpopulated throughout the Galactic Centre, which 
is not what is observed. We also attempted to reproduce 
the Galactic Centre cluster using an exponentially declining 
star formation rate, starting from a maximum at 14 Gyr. 
We found that, if the timescale on which the star formation 
rate drops off is significantly less than the assumed age of 
the cluster (14 Gyr), the model produces too few bright 
giants and too many faint ones. We concluded that, if the 
star formation rate has been declining, it must have done 
so slowly, so that the assumption of a constant rate is 
reasonable. 

In summary, the simple and conservative assumptions 
we have made about the IMF and star formation history 
of the Galactic Centre generate a model which matches the 
observed numbers of giants outside the depleted regions 
well and the alternative models we studied produced 
significantly worse fits. In Figure [2] we show the results of a 
Monte Carlo realisation of the Galactic Centre population 
constructed using these assumptions. The area of the 
square at each point in the (mass-K-magnitude) grid 
represents the relative number of stars of that mass which 
will be visible at th at K magnitude. W e see that the three 
bands observed by IGenzel et al.l l|l996l ) are dominated by 
low-mass objects in the range 1 — 4Mq. There is a small 
contribution from higher-mass objects, but it is suppressed 
both by the IMF, which makes these objects intrinsically 
rare, and by their rapid evolution. Middle band objects 
(12 < K < 10.5) are likely to be 1 - 2M Q stars, while 
bright band objects (K < 10.5) are likely to be 2 — 3Mq 
stars. Figures [T] and [5] clearly show that the missing objects 
are giant stars with masses in the range 1 — 3M© (the 
turnoff mass in our models is ~ 0.9 Mq). We now seek to ex- 
plain why there should be a deficit of these stars near Sgr A* . 



2.2 Possible causes of the depletion of giants 

It has been sug gested (e.g. iDavies fc Kin3l2005h that stars 
on the RGB or AGB on eccentric orbits about Sgr A* could 
be tidally stripped of their envelopes as they pass through 
periapse. Since objects on eccentric orbits spend most of 
their time at apapse, this could potentially produce a de- 
pletion of such stars at distances from Sgr A* much larger 
than the separations required to tidally strip RGB or AGB 
stars. We think that tidal stripping of this nature will make 
at most a minor contribution to the depletion of giant stars 
in the Galactic Centre, since stripping must occur during 
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the short RGB phase. This will form the subject of a later 
paper. 

Alternatively, numerous authors l|Genzel et al . 1996; 
iBailev fc Davieslll999l ; lAlexander|[l999l ') have suggested that 
the depletion of red giants is due to collisions of the giants 
with other members of the Galactic Centre stellar popu- 
lation. Collisions would remove mass from the giants' en- 
velopes and potentially alter their subsequent evolution so 
that they do not become bright enough to be visible in 
the depleted K bands. Since the Galactic Centre cluster is 
very densely-populated and has a high velocity dispersion 
(> 100 kms -1 ), collisions provide an attractive explanation 
for the lack of bright giants within the ce n tral ~ 0.1 pc. 

Previous work jDavies et all 1 19981 ; I Bailey fc Davies] 
1 19991 ; lAlexanderl 1 19991 ) has generally focussed on collisions 
impacting stars while they are on the RGB, and we largely 
agree that this is where the focus should be. It is highly 
likely that collisions of low mass MS stars with compact ob- 
jects could expel enough mass to take the stars below the 
turnoff mass (~ 0.9 Mq), ensuring that they never evolve to 
become giants within a Hubble time. However, this cannot 
be the sole cause of the depletion of giants in the Galac- 
tic Centre, as some giants are observed there. In addition, 
the effects of collisions on a dense cluster of main sequence 
stars are complex, since most such collisions are likely to be 
with other MS stars. Some fraction of these encounters will 
lead to mergers resulting in the conversion of two stars to 
a single higher-mass (and therefore brighter, but also more 
shortlived) star. This coagulation process fe.g. lLeell993h and 
the effect on MS stars of collisions with stellar-mass black 
holes will form the subject of a companion paper. 

As will be explained in more detail in a later section, 
we have modelled the effect of mass loss on stars at various 
different stages in their evolution using the STARS code. 
The mass losses required to substantially alter the evolution 
of HB or AGB stars are very large - > 99% of the enve- 
lope - and extremely unlikely even with a massive compact 
impactor. In addition, if a star is allowed to survive until 
the HB or AGB, it will already have been visible in the de- 
pleted K-bands for some of the duration of its RGB phase. 
Collisions occurring after the RGB will have very little effect 
on the observed stellar population; to alter the evolution of 
post-MS star so as to have a significant impact on the visi- 
ble population, collis i ons m ust occur on the RGB. 

IBailev fc Davies! (|1999T ) considered collisions between 
2 Mq and 8 Mq giants and impactors of a variety of masses 
intended to represent MS stars, neutron stars and white 
dwarfs. For each giant-impactor pair, they explored the 
(woo, Rmin) parameter space relevant for the Galactic Cen- 
tre. Throughout this paper, Voo denotes the relative velocity 
between the colliding objects at infinite separation and i? m i n 
is the minimum distance that would be reached by their cen- 
tres of mass, were they point masses on Keplerian hyperbolic 
orbits. 

Although IBailev fc Daviesl (|l999l ) did not quantify how 
much mass either of their two giants must lose to prevent 
them evolving to become visible in the middle or bright 
bands, they observed that the mass losses resulting from 
encounters of giants with MS stars, white dwarfs and neu- 
tron stars were not large, of the order of a few to ten percent. 
They concluded that such collisions were not able to explain 
the observed depletion of giants. 
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Figure 2. Plot depicting the relative contributions of stars of 
mass 1 — 12 Mq to K bands in the range 15 — 7, as seen at the 
Galactic Centre. The area of the squares is proportional to the 
number of stars of each mass in each K-band assuming a Miller- 
Scalo IMF, a constant star formation rate over 14 Gyr and using 
the Eggleton stellar evolution tracks and colours and bolometric 
corrections from Johnson 1966. The boundaries of the K-bands 
observed by Genzcl ct al. (1996) are indicated by the dashed hor- 
izontal lines. The large concentration of objects with K=15 — 16 
is the red clump. 



In this paper, we extend the work of IBailev fc Daviesl 

|l999l ) by considering, in addition to impacts with MS stars, 
encounters between giants and stellar-mass black holes. We 
also consider the effects of multiple impacts involving a 
single giant. We examine the effects of collisions involving 
1 — 3Mq giants and black holes or MS stars. 



3 COLLISION RATES IN THE GALACTIC 
CENTRE 

Within a few pc of the Galactic Centre, th e stellar dis- 
tribution follows a brok e n power-law cusp (|Genzel et al.l 
l2003l ; ISchodel et alj|2007jh iGenzel et all (|2003l ) find that the 
three-dimensional number density of visible stars n* varies 
with Galactocentric radius as 



-1.4 
-2.0 



for 
for 



r ^ 0.4pc 
r > 0.4pc 



n»(r) oc 
and lSchodel et al.l (|2007f > find 



for 
for 



r sC 0.22pc 
r > 0.22pc 



(1) 



(2) 



Since the break radii of both of these power-law models com- 
fortably exceed the radius within which depletion of bright 
giants is observed (~ 0.1 pc), we adopt a single power-law 
model for the visible stars of the form 



n,(r) oc 



(3) 
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Strictly adopting the power law of ISchodel et~ai1 (120071 ) af- 
fects only slightly the rates of collisions between giants and 
MS stars, and has little influence on our results. 

Normalising Equation [3] at Galactocentric radii of ~ 
0.1 pc is difficult, since the enclosed mass at this radius is 
dominated by the SMBH. The influence radius ri n fl, defined 
as the radius within which the enclosed stellar mass is equal 
to th e mass of Sgr A* is of the order of 1 pc fe.g lFreitag et all 

l2ooeh . 

To obtain a realistic normalisation at small radii, we use 
th e results of simulatio ns of the Galactic Centre performed 
bv lFreitag et all (|2006l ) using the ME(SSY)* *2 Monte Carlo 
stellar dynamics code l|Freitag fc B enz 2001, 2002;). The ini- 
tial conditions were r?-models l|Dehnenlll993l ) with rj = 1.5 
and a central massive object representing Sgr A* (mod- 
els with r\ = 2, corresponding to a shallower cusp, were 
also run, but produced very similar results). The mass of 
the central object, the mass of the Galactic Centre clus- 
ter and the break radius at which the 77-model steepens to 
p(R) oc R- 4 were taken to be 3.5 x 10 6 M Q , 7 x 10 7 M Q 
and 14 pc respectively , designed to match obs ervations of 
the G al actic Centre bylSchodel et al.l (|2003T ) and lGhez et al.l 
(|2005h . iFreitag et all | |2006l ) found that the MS stars relax 
into a cus p with 7 = 1. 3 — 1.4 , in agreement with the obser- 
vations of Genzel et al.l ((2003) . Using their results (e.g. their 
Figure 2), we adopt a normalisation for the main sequence 
stars such that M(O.lpc) = 4 x 10 4 M Q . For an evolved 
stellar population with a Miller-Scalo IMF, the mean stel- 
lar mass is ~ 0.5 Mq, so this normalisation corresponds to 
8 x 10 4 main-sequence stars within a Galactocentric radius 
of 0.1 pc. 

The population of stellar-mass black holes in the Galac- 
tic Centre is expected to have the form of a cusp with 
a power-law exponent o f 1.75 provided that they domi- 
nate the relaxation rate l|Bahcall fc WolJll976l ). This rela- 
tion sh ould hold at least in the in nermost few tenths of a 
parsec l|Amaro-Seoane et al.l [20071 . and references therein). 
However, setting the normalisation of this population is very 
diffic ult, since the black holes are not directly observable. 

iMorris! (|1993l ) estimated that all the black holes formed 
within ~ 4pc of the Galactic Centre would sink into the 
centre over the age of the Galaxy, resulting in a total mass 
in stellar black holes there of ~ 1O 6 M0. They were not 
able to say whether such a migration would result in a 
dense c luster of black holes, or would s imply add mass to 
Sgr A*. lMiralda-Escude fc Gouldl (|2000h estimated, in good 
agreement with IMorris! 1 19931 ). that all black holes formed 
within 5pc would migrate in over the age of the Bulge. Us- 
ing a model of the stellar mass densit y as a function of 
Galactocentric radius from lGenzel et all (|2000T ) and assum- 
ing a piecewise IMF from which all stars in the mass range 
30 - 100 M yield black holes, they estimated that 2.4 x 10 4 
BHs should have migrated into the central 0.7 pc. They cal- 
culated that the timescale on which stellar-mass BHs would 
be captured by Sgr A* would be considerably longer than a 
Hubble time, so that all these stellar-mass BHs should still 
be within this volume. A similar indirect means of estimat- 
ing the number of black holes within a given Galactocentric 
radius is given by lAlexander fc Livid <|2004h ■ They assume 
that stellar-mass BHs have migrated into the Galactic Cen- 
tre by mass segregation and estimate the maximum number 
that can survive there by requiring that the timescale on 



which the BH population is lost by scattering into the loss 
cone is longer than a Hubble time. For a population of 10 Mq 
BHs, their analysis yields a maximum number of BHs that 
can be within 0.1 pc of Sgr A* of ~ 5000. 

Stellar-mass black holes can be directly observed if they 
are accretin g material at suffici ent rates to be visible as X- 
ray sources. 



Muno et al l l|2005h reported four X-ray tran- 



sients within 1 pc of Sgr A*, an overabundance of a fa ctor 
~ 20 when compared to the field. lBaganoff~e t al. ( 2003|) de- 
tected diffuse X-ray emission with a 2 — 10 keV luminosity of 
10 33 ergs -1 coincident with the position of Sgr A* to within 
0.27arcsec and with an intrinsic size of 0.61 arcsec. They 
consider the possibility that the source may be a population 
of compact objects accreting from the ambient gas, but point 
out that the luminosity of such a system would be orders of 
magnitude less than that of Sgr A* itself accreting from 
the same gas. They point out that a single X-ray binary, in 
which the compact object is supplied with a strong accretion 
flow, would have a luminosity comparable to that observed. 
However, the number of X-ray binaries expected within this 
volume (r < 0.03pc) is very small, since binaries are diffi- 
cult to form in a region with such a high velocity dispersion, 
and would also be quickly disrupted by stellar encounters. 
Taking into account that the observed extent of the source is 
consistent with Sgr A*'s Bondi radius, ffBaganoff et al.l (|2003l ) 
conclude the emission is most likely to be due to accretion 
ont o the supermassive bla c k hol e itself. 

iDeegan fc Navakshinl l|2007h examined the question of 
whether the X-ray emission could be due to compact ob- 
j ects accreting from the ISM, specifically the Minispiral 
jPaumard et alll2004 ). in more detail. They modelled the 
Minispiral as a half-disc of 50 Mq of gas in a circular orbit 
between 0.1 and 0.5 pc around Sgr A*. They constructed a 
Bahcall and Wolf cusp of black holes extending to 0.7 pc, al- 
lowed the black holes to accrete the gas and calculated their 
instantaneous luminosity as they moved along their orbital 
tracks. By following how many black holes were detectable 
at any given time, they placed constraints on the maximum 
number of compact objects. They found that the number 
of stellar-mass BH within 0.1 pc was likely to be at most a 
few thousand. None of these authors consider the possibil- 
ity that X-ray sources may be lone stellar-mass black holes 
accreting material captured during a stellar colli si on. 

We again use the results of IFreitag et al.l (|2006l ) to 
obtain a fiducial model for the number and distribution 
of black holes in t he innermost few tenths of a parsec. 
IFreitag eT al. (2006) found that the stellar mass black holes 
migrated inwards by mass segregation, for ming a cusp with 
7 = 1.8, in very good agreement with iBahcall fc WolJ 
(|l976l ). and that the central 0.1 pc contained a fewx 10 3 BH. 
We therefore adopt a black hole population with a cusp of 
the form 



n B H(r) 



(4) 



normalised so that riBH(O.lpc) = 5 x 10 3 . The radius at 
which the total enclosed mass of our MS and BH populations 
becomes equal to the mass of Sgr A* (4 x 10 6 M Q ) is then 
~ 1.4pc 

In this work, it will be assumed that stars lie inside the 
sphere of influence and hence that their orbits are Keplerian 
with one-dimensional velocity dispersion <tid is oc r -1 / 2 . If 
the stellar density distribution has the form of a cusp with 
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three dimensional radial distribution oc r 
dimensional velocity dispersion is given by 



o"id(f) 



GM 



(7 + l)r 



then the one- 



(5) 



1000 



If a star of mass M» and radius 7i* moves through a 
cluster of impactors of mass Mi mp and radius i?i mp such 
that the relative encounter velocity is Voo , the effective cross 
section s of the star is given by 



s = ir(R* + Riv 



2G(M, + M im( 
(R* + R imp )vi 



(6) 



The second term inside the braces on the right-hand side of 
Equation[6]is the increase in the star's effective cross section 
due to gravitational focussing. 

At a given Galactocentric radius r, stars will in reality 
have a distribution of velocities which must be integrated 
over to yield an accurate determination of the collision rate 
at that radius. We find that the velocity distribution in an 
?7- model is very similar t o a Ma xwellian distribution. Follow- 
ing [Binne^&_I¥emaii3 (|l987h . we define the collision time 
averaged over the velocity distribution by 

icoiM" 1 = 4 v / 7r(ii» + i?imp) 2 n imp (r)CTiD(r)x 



1 + 



G(M»+M in 
2a'f- D (R,+R i 



np)J 



(7) 



where ni mp (r) is the number density of the impactors. 

Since the radius of a star changes drastically during its 
lifetime, R, — R*(t) and t co \ is a function of time as well as 
Galactocentric radius. As a result, a more useful quantity 
is the integrated probability p(r) of the star experiencing a 
collision over its lifetime, given by 



p(r) 



dt 



t C oi(r,t) 



(8) 



where 7* is the lifetime of the star. In this paper, we are pri- 
marily interested in collisions occurring to stars while they 
are on the giant branch, in which case the collision proba- 
bility of interest is 



HCTRGB) dt 
PRGB(r) = / 

J £(BRG 



j t col (r,t) 



(9) 



where r,(BRGB) is the age of the star at the base of the 
red-giant branch and i(TRGB) is its age at the tip. We 
define the base of the red-giant branch as the point at 
which the star's luminosity reaches a minimum at the 
end of the Hertzprung gap, and the tip of the red-giant 
branch as the point where the star's luminosity reaches its 
maximum value on the giant branch. 

In Figure |31 we plot the integrated collision probability 
given by Equation [9] as a function of Galactocentric radius 
for 1, 1.5 2 and 3Mq stars, assuming that the impactors 
are 1 Mq main-sequence stars (solid lines) or 10 Mq black 
holes, distributed according to Equations [3] and [4] For 
reasons of clarity, we have assumed that the inner power 
law a = 1.4 holds at a ll radii - sinc e the b reak in the cusp 
power law derived by iGenzel et al.l ()2003h occurs at radii 
where the integrated probability of suffering a collision with 
a MS impactor is small (< 10%), this does not affect our 
results. 

In deriving Equation |9j it is assumed that ni mp (r) 



a 10 ; 



0.1 ET 



0.001 = 



_ _ RGB-BH 



0.001 0.01 0.1 1 

Galactocentric radius (pc) 

Figure 3. Plot showing the integrated probability that stars 
of masses (from top to bottom) 1, 1.5, 2 and 3Mq will collide 
with main-sequence stars (solid lines) or stellar-mass black holes 
(dashed lines) while themselves on the red giant branch, as a 
function of Galactocentric radius. 



and Voo axe constant over the lifetime of the star, which 
can only be true if either the star's orbit is circular, or if 
the star's lifetime is much shorter than its orbital period. 
Since neither of these of these conditions is true in general, 
we integrated the collision probability self-consistently on 
orbits with a variety of eccentricities. The effect of orbital 
eccentricity was to increase the probability of a star with 
a given semimajor axis suffering a collision by factor of 
at most 2 over the collision probability obtaining from a 
circular orbit with the same semimajor axis. 

A further simplification made in estimating these 
collision rates is that the impactors are all either 1 Mq 
main-sequence stars or 10 Mq black holes. Given the lack 
of knowledge about the mass function of stellar-mass black 
holes, the latter assumption is reasonable (and will be 
discussed in more detail later), but the former may not be. 
We therefore recalculated the collision rates assuming that 
the main-sequence impactors were distributed according to 
a Miller-Scalo IMF between 0.08 and 1M . We found that 
this makes little difference to the collision rates. 

We see from Figure[3]that the probability of a 1 — 2 Mq 
giant suffering a collision with a main-sequence star exceeds 
unity within a radius of ~ 0.1 pc and the probabilities of 
the same giants being struck by a black hole exceeds unity 
within a radius of ~ 0.03 pc. The 3Mq star has a very 
short giant phase in comparison to the lower-mass objects, 
and the radii at which the probabilities of the 3 Mq giant 
being struck approach unity are approximately an order of 
magnitude smaller. 

We interpret collision probabilities in excess of unity 
as implying that stars at that radius will suffer more than 
one collision. Figure [3] shows that the region of space within 
0.1 pc of Sgr A* is a highly collisional environment. In 
the following sections, we investigate the effect of stellar 
collisions on the population of bright giant stars in the 
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Galactic Centre. 



particle grid was then deformed according to 



rt = 



4 NUMERICAL METHODS 
4.1 Stellar evolution code 

To evolve the stellar models presented in this paper we used 
STARS, the Cambridge St ellar Evolution code. STARS was 
originally written by Peter lEggletonl (1 19Tlh an d has been ex- 
tensively modified since fsee lPols et al.lll995l and references 
therein for a complete de scription). We use the Re imers 
mass-loss law o n the RGB llKudritzki fc Reimersl[l978l ) with 
r\ = 0.4 and the lVassiliadis fc Wood! 1 19931 ) mass-loss law on 
the AGB; on the MS the mass-loss is assumed to be neg- 
ligible for the low-mass stars considered h ere. Convective 
overs hooting was included with 8 ov = 0.12 (Sc hroder et al.l 
1 19971 ). 

Following lEldridge fc Vmk] ( 2006) we modify our mass- 
loss rates for high -mass stars. For O B stars, we use the mass- 
loss predictions o f lVink et all (|200lh which scale with metal- 



licity as (Z/Zc. 



vl/2 



For all other pre -WR phases, we em- 



ploy the rates of Ide Jager et al] (|l988l ) scaled similarly with 
metallicity. When the star becomes a WR star (A sur f < 0.4 
and lo g(T/K) > 4.0), we use the rates of iNugis fc Lamersl 
i|2000l ). 



4.2 Smoothed particle hydrodynamics 

For the hydrodynamical simulations presented here, we 
used a Smoothed Particle Hydro dynamics (SPH) code 
based on that described in iBenzl (jl990l ). The code uses 
the SPH formalism to solve fluid equations and a binary 
tree to calculate gravitational forces. The code has been 
modified to include point masses, particles which only 
interact with other particles via grav itational forces an d 
are used to model compact objects l|Davies et al.l 1 199 ll ). 
Their gravitational fields are smoothed in the same way as 
those of ordinary SPH particles, although the smoothing 
lengths of the point masses are fixed. We use the standard 
SPH artificial viscosity formalism with a = 1, /3 = 2. Gas 
particles are assumed to behave adiabatically. 

We constructed SPH stellar models using one- 
dimensional density and temperature profiles from the 
STARS code. SPH particles are placed on a uniform 
hexagonal close-packed grid and their masses iteratively 
adjusted so that the local density matched that of the 
one-dimensional model. 

The extreme density gradients in the envelopes of 
giant stars, both very near the core and near the stellar 
surface, are difficult to model. In order to reproduce the 
density profiles of the outer envelopes of our giants while 
retaining sufficient resolution of the material near the core 
and avoiding very large particle numbers, we deformed the 
uniform particle grid in our models. Using results from our 
stellar evolution calculations, we determined the critical 
mass loss required to prevent the giant evolving to become 
bright enough to appear in the depleted K-bands. This 
gave us a minimum mass (enclosed by a radius r cri t) that 
we needed to resolve. We constructed a uniform particle 
grid so that r cr u contained at least 2000 particles. The 



ricosh a 



I j | for n > r clit 
for n < r crit 



(10) 



where n is the initial radius of a given particle, rf is the 
final radius and a is a constant ensuring the outer radius of 
the deformed grid is equal to the radius of the star. 

We modelled the giant cores as point masses. To 
stabilise the interaction between the point mass core and 
the envelope, the innermost 477 gas particles (i.e. those 
surrounding the core, corresponding to an enclosed mass 
about 1% of the envelope) were 'frozen' - their velocities 
were constrained during simulations to be equal to that of 
the core particle. Ensuring adequate resolution of the outer 
envelope and the core region and guaranteeing that the 
model be stable required ~ 175, 000 particles. 

The models were allowed to relax in isolation for 
several freefall times with a linear damping term acting 
to remove any residual oscillations in the particle grids. 
Collision simulations were started with the stars sufficiently 
far apart that tidal effects were negligible. 

Most calculations were run on single cpus on PCs. 
However, in order to examine numerical convergence, we 
repeated several simulations at high-resolution on the 
UK Astrophysical Fluids Facility at Leicester Univer- 
sity, UK. In these calculations, we used ~ 1.2 million 
particles in our giant models. The factor of ~ 7 greater 
particle number gives a factor of ~ 2 better linear resolution. 



4.3 Restricted three— body code 

In order to extend the range of simulations we could perform 
without consuming unacceptable quantities of computer 
time, and as an independent check on our SPH calculations, 
we also studied red giant collisions using a simple restricted 
three body code. As pointed out by iLivne fc TuchmarJ 
( 1988), if the impact velocity is much greater than the speed 
of sound in the giant's envelope, the envelope cannot react 
to the passage of the intruder. The interaction can then be 
approximated by an encounter between two point masses 
(the impactor and the giant core) with the envelope repre- 
sented by a third, extended particle with the enclosed mass 
profile of the undisturbed envelope. The force between the 
core and the intruder -Fc-int is then given by 



GM COIC M h 



(11) 



and the force between a point particle i and the envelope 
particle is given by 



Fi- 



GMiMenvCn- 



(12) 



where M cnv (ri_ cnv ) is the envelope mass enclosed within a 
radius ri_ env . The equations of motion of the three particles 
can then be easily integrated. 

In collisions in which the impactor imparts a significant 
velocity to the core relative to the envelope, the fraction 
of the envelope retained by the core can be estimated by 
computing the mass within the Bondi-Hoyle radius, given 
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by 



GM a 



(13) 



where v TO \ is the relative velocity between the core and the 
envelope. This analysis can be improved upon slightly by 
instead computing the mass contained within the Bondi- 
Hoyle-Lyttleton radius, given by 

GMcorc 



»~BHL = 



+ Cs(r B HL) 



(14) 



where tbhl must be found iteratively since it depends on 
the sound speed at that radius. We find in practice that the 
two schemes give very similar results. 

The envelope mass retained by the core may thus be 
estimated for any given collision by using the three-body 
technique to find the relative velocity induced between 
the core and the envelope, yielding rBH or tbhl, and 
determining the enclosed mass at this radius from the 
giant's density profile. 



4.4 Monte Carlo multiple— collision code 

Since the collision rates shown in Figure [3] imply that 
giants may suffer multiple collisions with a mixture of 
black-hole and main-sequence impactors, we explored this 
possibility using a Monte Carlo technique. For a population 
of 10 4 giants, we simulated collisions at all the relative 
velocities for which we have SPH simulations (400, 800 and 
1200kms~ 1 for the 1 and 2M giants, 800kms _1 only 
for the 1.5 and 3Mq objects) and at randomly-chosen 
periastrons (distributed so that the probability of a collision 
at a periastron 7? is proportional to R 2 , appropriate for 
non-focussed encounters) with the numbers of black hole 
and main-sequence impactors being set by the collision rate, 
and thus the Galactocentric radius. For each impact, we first 
constructed a synthetic mass-loss-against-7? m i n curve from 
the results of the SPH and three-body simulations of that 
encounter. As we will show in detail later, the three-body 
code agrees very well with the SPH calculations in encoun- 
ters at small periastrons, in which the core is ejected from 
the giant, but hydrodynamic effects become important at 
larger periastrons (particularly at lower relative velocities). 
We therefore use the smooth three-body curve for very 
close encounters where our SPH calculations are sparse, 
out to a periastron where it begins to underestimate the 
mass loss derived from the SPH calculations, and from this 
periastron outwards, we use the SPH mass-loss curve. Using 
this synthetic mass-loss curve, we calculate the mass lost 
at the randomly-chosen periastron by interpolation. The 
mass lost in successive collisions can either be assumed to 
be multiplicative or additive. If a given collision removes a 
fraction of mass /i, the total mass lost f to t after N collisions 
can be estimated by 1 - / tot = (1 - /i)(l - /a)-»(l - /at), 
or by /tot = fi + ,f'2 + ---+fN- We found that it made 
little difference to the results which of these assumptions 
was made. From the total mass lost by each giant, we 
determined the fraction of the 10 4 objects prevented from 
evolving into the bright and middle bands for each impactor 
number. 



5 MODELLING THE EFFECT OF 

COLLISIONAL MASS LOSS ON THE 
EVOLUTION OF STARS 

To study the effects of mass loss on the giants, we con- 
structed in the STARS code models of stars of 1, 1.5, 2 
and 3Mq at several different ages on the main sequence, 
red giant branch, horizontal branch and AGB. 

In order to simulate the rapid mass loss caused by a 
collision we removed the desired quantity of mass from the 
star at a rate of 5 x lO _8 M0yr _1 , which is the largest mass- 
loss rate at which we could reliably converge models for 
low-mass stars. This rate is much lower than the mass-loss 
rates observed in our SPH calculations, which could reach 
~ lO 4 M yr _1 (the equivalent of ~ lM hr _1 )but is still 
much faster than the evolutionary timescales for low-mass 
stars. During this time we turned off composition changes 
owing to nuclear burning but retained the energy generation; 
this prevents nuclear evolution occurring whilst the mass is 
removed but also keeps the star's structure consistent. This 
is equivalent to making the collisional mass-loss instanta- 
neous on a nuclear timescale. 

We removed increasing fractions of the mo dels' en- 
velope s and studied their evolution, using the Ijohnsor] 
(1966) colours and bolometric corrections to calculate the 
K-band magnitude of the objects as functions of time. We 
were therefore able to determine, for a star of a given mass, 
how much mass a collision at any given stage in the star's life 
would need to expel such that the star would never evolve 
to become brighter than a given K-magnitude. In Figure 3J 
we show by means of four HR diagrams the effect on the 
evolution of a IMq star of removing increasing fractions of 
its envelope at a point halfway up the giant branch. In each 
panel, the normal evolution of the star is denoted by the 
thick black line. The star is then taken to instantaneously 
lose mass, leaving it with the envelope fraction given in each 
figure panel. We see that leaving the object with 93% of its 
envelope has little effect on its evolution - it follows a normal 
AGB track and becomes a white dwarf. However, leaving the 
star with 79% of its envelope results in an object resembling 
a subdwarf-B star, with a short AGB significantly fainter 
than that of the undisturbed star (although the luminosity 
of the tip of the RGB is not significantly changed). Leaving 
the star with 66% or 56% of its envelope, as shown in the 
lower panels of Figure [3] prevents the star igniting helium, 
thus removing the AGB phase altogether, and also makes 
the tip of the RGB fainter. In summary, removing a few 
tens of percent of the envelope of a low mass giant is suf- 
ficient to make its AGB phase significantly fainter, but to 
radically alter its evolution such that it never ignites helium 
at all requires the removal of ~ 40% of the envelope. 

In all post-main-sequence phases, the fractions of the 
stars' envelopes that must be expelled in order to signifi- 
cantly affect their evolution are substantial, and especially 
so for HB or AGB stars. We find that, in order to signifi- 
cantly alter the evolution of a star, mass must be removed 
whilst it is on the main-sequence or the red giant branch. 
In Figure [5] we illustrate this point by comparing the effect 
of removing 60% of a IMq star's envelope at three different 
phases of its life. Each panel in the plot shows the evolu- 
tion with time of the star's K-band magnitude. The top left 
panel shows the evolution of the undisturbed object. The 
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Figure 4. Plot showing the effect on the evolution of a IMq star of removing increasing fractions of its envelope mass, once it has 
reached a point halfway up the RGB. The star's normal evolution up to this point is shown by the thick black line, and the evolution 
after mass loss by the dashed line. The fraction of the object's envelope remaining is given in the bottom right of each panel. 



top-right panel shows the result of removing 60% of the en- 
velope mass halfway along the RGB, the bottom left panel 
shows the result of removing this mass halfway along the 
horizontal branch and the bottom right panel shows the re- 
sult of removing mass on the early AGB. In all cases, the 
solid lines represent evolution of the undisturbed object and 
dashed lines represent evolution after mass loss has occurred. 
The plots clearly show that this mass loss occurring on the 
RGB has a strong effect on the star's evolution, making the 
tip of the giant branch ~ 1.5 magnitudes fainter (although 
having little effect on the duration of the RGB phase), and 



preventing the object evolving into the HB or AGB phases. 
Conversely, removing the same fraction of envelope mass on 
the HB or AGB has almost no effect on the star's subsequent 
evolution and, in particular, does not noticeably reduce the 
brightness of the tip of the AGB. The effects of mass loss on 
main-sequence stars will be discussed in a companion paper 
and here we concentrate on mass loss due to collisions on 
the red giant branch. 

Stellar evolution on the RGB is very fast and stellar 
properties, particularly core mass and radius (and hence lu- 
minosity) change significantly on the RGB. The age of a star 
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Figure 5. Plot showing the effect on the evolution of the K— band brightness of a 1 Mq star of removing 60% of its envelope mass halfway 
along the RGB (top right panel), halfway along the HB (bottom left panel) and on the early AGB (bottom right panel). Top left panel 
and solid lines in other plots represent undisturbed evolution. 



is not a useful coordinate against which to measure its evolu- 
tion and instead, we define a quantity based on the evolution 
of the core mass, which we use as a proxy for time: 

( f M COTe (t) - M cole (BRGB) 

Tc ° rcU Moore (TRGB) — Mcore(BRGB) ( ' 

where M CO rc(BRGB) is the core mass at the base of the 
red giant branch and M core (TRGB) is the core mass at 
the tip. To illustrate the relationship between T core and 
the actual age of the giant, we plot r cor c as a function of 
time on the RGB for a 1 Mq giant in Figure [6] We show 
in Figure [7] the lifetimes of 1 Mq giant stars in the middle 
band (12 > K > 10.5) after encounters at five different 
positions along the giant branch (characterised by r CO rc on 



the y-axis) which leave them with between 1 and 100% of 
their envelopes remaining. 

The mass loss required to prevent the 1 Mq giant 
evolving to become visible in the middle band increases 
as the object ascends the giant branch. As the core mass 
grows, its later evolution becomes less and less sensitive 
to the envelope mass. However, even near the base of the 
giant branch, the quantities of mass that must be removed 
to seriously affect the giant's evolution are substantial - 
several tens of percent at least. 
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M(M ) 


M core (M Q ) 


R(R ) 


Age(yr) 


"7~corc 


K(mag) 


1.0 


0.30 


25.2 


1.22 x 10 10 


0.51 


14.5 


1.5 


0.29 


18.8 


2.95 x 10 9 


0.62 


15.5 


2.0 


0.34 


31.9 


1.21 x 10 9 


0.54 


13.4 


3.0 


0.39 


27.6 


3.76 x 10 s 


0.40 


13.6 



Table 1. Properties of giants used in our SPH calculations. 
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Figure 6. Plot showing the relationship between t co 
elapsed time on the RGB for a 1 Mq giant. 



and 



6 COLLISIONS INVOLVING GIANT STARS 



6.1 



Collisions of giants with stellar— mass black 
holes 



The depleted K-bands reported bv lGenzel et alj (|l996t ) are 
dominated by 1, 1.5, 2 and 3 Mq stars, as shown in Figure[3] 
We chose to model collisions involving giants approximately 
halfway up the giant branch, as measured by the quantity 
T corc , taking these objects to be representative of giants of 
a given mass. This has the advantage that such objects are 
not yet bright enough to be visible in the depleted K-bands 
and that they lie close to the point at which the cumulative 
probability of suffering a collision while on the giant branch 
is approximately 0.5. We give the parameters of our model 
giants in Table 1. 

The quantities of interest in these simulations are the 
masses of gas bound to the giant core, and captured by 
the intruding black hole. To calculate these masses, we 
employed an iterative scheme. For each point mass, the 
gas particles in its immediate vicinity were fetched, their 
velocities relative to the point mass determined and hence 
their kinetic, potential and thermal energies in the frame 
of the point mass were calculated. Those particles with 
negative total energy were taken to be bound to the point 
mass. A composite object was then constructed from the 
point mass and those particles bound to it. The mass, 
centre of mass and centre-of-mass velocity of the composite 
object were calculated and the procedure repeated until the 
mass of the composite object converged. 




Fraction of envelope mass remaining 



Figure 7. Plot showing the effect on the timescale for which a 
IMq giant is visible in the middle band (12 >K> 10.5) of strip- 
ping mass from the giant at various stages of evolution on the 
giant branch. The evolutionary stage at which mass is stripped 
is denoted by the increase in core mass as a fraction of the total 
change in the core mass along the whole giant branch (defined in 
Equation 1151) , given on the y— axis and the fraction of the enve- 
lope mass remaining after stripping is given on the x— axis. The 
symbols denote the length of time for which each stripped model 
is visible in the middle band for the duration of its life. 



After a sharp decrease around the time of the black 
hole impact, the mass retained by the giant core declined at 
an ever-decreasing rate for long periods of time. To avoid 
running simulations for prohibitively long times, we devised 
a procedure to estimate the final mass of the collision 
product. Plotting the fraction of the envelope bound to the 
core, / cnv as a function of 1/t revealed that, except near 
the time of the impact, / env could be well fit by a straight 
line and extrapolated by least squares fitting to 1/t = 
(i.e. t — oo), yielding the final envelope mass retained by 
the core. In Figure [H] we show an example plot, taken from 
the collision of a IMq with a lOM© black hole at an J? m i n 
of IORq and of 800km s . This plot shows that at 

the late stages of the collision (i.e. for small 1/t), the curve 
is very well fit by a straight line, giving a robust estimate 
of the final mass loss from the encounter. 

To check for numerical convergence, we repeated several 
calculations at much higher resolution on the UKAFF facil- 
ity at the University of Leicester. In Figure [9] we compare 
the evolution of the orbital separation (top panel) of the 
giant core and impacting black hole, and of the energies of 
the system (lower panel), in the low-resolution calculation 
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Figure 8. Plot of the envelope fraction remaining against 1/timc 
for the an encounter between a 1 Mq giant and a 10 Mq black 
hole with a v^x, of 800 kms -1 at an -R m i n of 10 Rq. The solid 
line with the crosses is the measured envelope fraction remaining 
from the SPH calculation and the dashed line is the fit applied to 
estimate the final remaining envelope fraction. 

(black squares) and the high-resolution calculation (solid 
lines) of an encounter between a 1 Mq giant and a 10 Mq 
black hole with a v x of 400 km s -1 at an i? m i n of 10 Rq. 
The agreement between the two calculations on the gross 
evolution of the system is clearly excellent. 

In Figure [TTJ1 we show density slices in the 2 = plane 
1.6 crossing times after periastron passage. The agreement 
between the low-resolution (left panel) and high-resolution 
(right panel) resolution calculations is again very good. The 
same structures are visible in both images, including the 
disk of captured material around the black hole, and the 
wake of the black hole through the red giant, including the 
hole where it has exited the envelope. We also find that the 
low- and high-resolution runs agree well on the fraction 
of the envelope retained by the core, giving 15% and 17% 
respectively in these calculations. 

We observed, particularly at smaller periastrons, that 
the mass loss was often due to the core being ejected and 
carrying part of the envelope away with it. This is clearly 
vis ible in Figure 1101 This phenomenon was pointed out 
by iLivne fc Tuchmanl (|l98Sf ) and is essentially due to the 
inability of the envelope to react to the supersonic impact of 
the intruding black hole. We use our restricted three-body 
code to determine the velocity of the core relative to the 
envelope and estimate the mass retained by the mass within 
the Bondi-Hoyle-Lyttleton radius corresponding to that 
velocity. 

For each of our SPH giant-black hole encounters at 
a given v-a, we performed reduced three-body calculations 
over a range of periastrons and calculated the mass retained 
by the core using the Bondi-Hoyle/Bondi-Hoyle-Lyttleton 
formalism. In Figure 1111 we show the results of SPH 
calculations of collisions between a 1 Mq giant star and a 
10 Mq black hole performed at a range of periastrons and at 
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Figure 9. Plot comparing the evolution of the separation between 
the red giant core and the black hole (top panel) and the energies 
of the system (bottom panel, from top to bottom, kinetic, ther- 
mal, total ad gravitational potential energies) of high-resolution 
and low— resolution simulations of a collision between a 1 Mq gi- 
ant star and a 10 Mq black hole at a periastron of 10 Rq and 
Voo of 400kms _1 . Solid lines are results from the high— resolution 
calculation and black squares are from the low-resolution calcu- 
lation. Note that the initial separation between the giant and the 
black hole in the high-resolution calculation was smaller to save 
computer time. 



Woo = 400, 800 and 1200 kms -1 (black lines with squares) 
compared with the results of our reduced three-body 
simulations over a similar range of periastrons and with 
Voo = 800 km s" 1 . 

The SPH simulations and the reduced three-body 
analysis agree very well at small periastrons (< 10 Rq), 
but depart at larger values, with the SPH results giving 
systematically smaller fractions of the envelope being 
retained (i.e. predicting that more mass is expelled). 

At the high velocities considered in this paper, there 
are two regimes of mass loss in interactions between giants 
and stellar-mass black holes, delineated by the periastron 
of the collision. At large periastrons (> i? g i an t/2), mass 
loss is largely due to the shock driven by the impactor 
expelling parts of the giant's outer envelope, leaving the 
core region relatively undisturbed. Conversely, at perias- 
trons < -Rgiant/2, the core receives a strong impulse and is 
displaced from the envelope, carrying some of the envelope 
material away and leaving the rest of the giant to disperse 
on a dynamical timescale. 

The critical periastron, inside which collisions expel 
enough mass to prevent the 1 Mq giant evolving to become 
brighter than K — 12 at the Galactic Centre, is predicted 
to be ~ 12 Rq by the SPH calculations and ~ 10 Rq 
by the Bondi-Hoyle analysis, resulting in a significant 
difference in the critical collision cross section if collisions 
are unfocussed. Hydrodynamic effects thus cannot be 
neglected in this problem and SPH simulations are essential 
in determining the collision parameters required to prevent 
giants evolving into the brighter K-bands in the Galactic 
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Centre. 

Since little is known about the mass function of black 
holes, we explored the consequences of impacts with black 
holes of different mass, to see if, for a black hole population 
of a fixed total mass a smaller number of higher-mass black 
holes would have a greater affect on a population 1 Mq 
giants than a larger number of lower-mass black holes. For 
this to be the case, the critical encounter periastron Rcriticai 
for stripping enough of the giant's envelope to prevent it 
becoming visible in the depleted bands must increase faster 
than -Rcriticai oc \^Mbh, if collisions are not focussed. We 
found that this is not the case and hence that the increase 
in -Rcriticai with black hole mass was not enough to offset 
the smaller number of impactors and therefore that a 
population of black holes with the canonical mass of 10 Mq 
is likely to have the greatest effect on the giant population. 

In collisions with 10 Mq black hole, the critical pe- 
riastron required to prevent the 1 Mq giant evolving to 
become brighter than K — 12 is ~ -Rgiant/2. Hence, if 
collisions are unfocussed, ~ 25% of single encounters will 
prevent the giant evolving to become brighter than K — 12 
at the Galactic Centre. This fraction is, however, a lower 
limit. The collision rates depicted in Figure [3] imply that 
giant stars can expect to experience > 1 encounter with 
stellar-mass black holes in the innermost O.lpc of the 
Galactic Centre, roughly the region where the population 
of giants is observed to be depleted. In addition, Figure [TT] 
shows that significant mass loss occurs at periastrons larger 
than the critical one, so several successive encounters may 
remove a large enough fraction of the a giant's envelope to 
affect it to the degree required. 

We used our Monte Carlo code, with black-hole im- 
pactors alone, to generate the fraction of giants prevented 
from evolving into the middle and bright bands for a given 
number of impactors at three velocities for the 1 and 2 Mq 
giants and a single velocity for the 1.5 and 3Mq giants. 
The results of our calculations for the 1 Mq giant are shown 
in Figure [T2l 

A given Galactocentric radius uniquely determines a 
relative velocity and cumulative impact probability. To map 
the fractions of objects prevented from evolving into the 
bright and middle bands, we performed two-dimensional 
interpolation in Voo and collision probability for the 1 and 
2Mq giants to determine the fraction of these objects 
prevented from evolving into the middle and bright bands 
as a function of Galactocentric radius. For the 1.5 and 3Mq 
giants, we interpolated in collision probability only. 

This analysis assumes that all giants of a given mass 
may be represented by the models at the halfway stage on 
the giant branch used in the SPH calculations. To check 
the influence of the evolution of the giants, we repeated 
the Monte Carlo simulations using 10 4 1 Mq giants with 
ages randomly distributed on the giant branch. Once the 
age of each giant was chosen, a three-body mass-loss/R m in 
curve was generated by interpolating between the mass-loss 
curves generated from the three-body calculations using 
giants bracketing the chosen age. As Figure [11] shows, the 
three-body calculations underestimate the mass loss at 
large values of i? m i n where hydrodynamic effects become 
important. We allowed for this by calculating from Figure 
[11] the factor by which the R m in resulting in a given fraction 
of the giant envelope was larger in the SPH results than in 



the three-body results. Assuming that these factors were 
the same for giants of any age, we produced a synthetic 
mass-loss curve for each giant of the randomly-chosen age. 
Collisions with these giants and a range of numbers of 
black hole impactors were then performed in the same way 
and the fraction of giants prevented from evolving into the 
middle and bright bands as a function of Galactocentric 
radius determined. In Figure [T3l we compare the fraction of 
1 Mq giants prevented from evolving into the middle band 
as a function of Galactocentric radius, taking all giants to 
be the representative T CO rc = 0.51 model (solid line), or 
accounting for the giant's evolution in the way described. 
We see that the giant's evolution has little influence on 
the results and that the models we chose to perform 
hydrodynamic simulations with are well representative of 
giants of each given mass. 

We see that collisions with black holes alone reduce 
the population of giants in the middle band by > 50% 
inside a radius of 0.04 pc, comparable to the size of the 
region depleted in this band. However, we will not draw 
any conclusions about the overall effect of collisions on the 
giant population until we have considered collisions with 
main sequence stars, described in the next section. 

In all of our SPH calculations, some material was 
captured by the black hole during its passage though the 
giant's envelope, the captured gas remaining in orbit in 
a disc-like structure around the hole. The black hole will 
accrete this material and the resulting energy release will 
make it visible as an X-ray source. In order to see how many 
such sources we might expect to observe, we assume that 
the captured material will be accreted at the Eddington 
rate, given by M = 2 x 10 7 M yr -1 for a 10 M black hole. 
The black holes typically capture a fewxlO -2 — 1O _1 M 
of material, so the time for which they will be visible as 
accreting X-ray sources is ~ 5 x 10 4 — 5 x 10 5 yr. If we 
assume that the collision rate per giant is given simply 
by riBHOgiantf co and take fiducial values appropriate for 
the innermost O.lpc of jibh = 10 6 pc -3 , a g iant = TrRgiant 
(neglecting gravitational focussing) with R g i a nt = 50 Rq 
and Voo = 1000 km s" 1 , we obtain a collision rate per giant 
of 4 x 10 -9 yr _1 . Our Galactic Centre models imply that 
the central O.lpc should contain ~ 10 4 giants, so that the 
rate of giant-BH collisions within this volume would be 
4 x 10 -5 yr -1 . Since this figure is approximately the inverse 
of the accretion timescale, we would expect to observe 
one or zero black holes accreting material which they have 
acquired durin g a collision with a giant. This result is 
consistent with lMuno et all l|2005h . 



6.2 Collisions of giants with main— sequence stars 

Encounters between giants and MS stars have already been 
studied in the context of th e giant depletion in the Galactic 
Centre by iBailey fc Daviesl (Il999l ) . How ever, the two giant 
stars studied bv lBailev fc Daviesl (1 19991 ) were both in quite 
extreme evolutionary stages. Their 2 Mq giant was near the 
tip of the RGB and their 8Mq giant near the tip of the 
AGB. We have chosen to study encounters with giants to- 
wards the middle of the RGB which are representative of all 
giants of a given mass, so we extend our study to include en- 
counters between giants and 1 Mq MS impactors, which we 
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(a) Low-resolution run (1.75 X 10 5 particles) (b) High-resolution run (1.19 X 10 6 particles) 



Figure 10. Comparison of density slices taken through the 2 = plane of snapshots from the low- and high— resolution runs 1.6 crossing 
times after periastron passage in the encounter between a 1 Mq giant star and a 10 Mq black hole at a periastron of 10 Rq and Voo 
of 400 kms -1 . Yellow colours represent densities of 6 X 10 — 2 gcm~ 3 and red represents densities of 6 X 10~ 5 gcm — 3 . The black hole 
(moving from right to left) and and the giant core are represented by white dots. Images are 150Rq on a side. The ejection of the core 
and retention of some of the envelope is clearly visible. In this calculation, the envelope fraction retained by the core was ~ 15% in the 
low— resolution run and ~ 17% in the high-resolution run, which is sufficiently small that the giant will not evolve to become visible in 
the middle band. 



a gain treat as point masse s. In addition, we extend the work 
of iBailev fc Daviesl (Il999l ) by considering the cumulative ef- 
fect of multiple impacts. We do not consider white dwarf 
or neutron star impactors. This is partly because, from a 
numerical point of view, there is little difference between a 
white dwarf, a neutron star and a low-mass main-sequence 
stars, since they are all point masses of similar mass. In ad- 
dition, our population synthesis models suggest that that 
MS stars outnumber white dwarfs by ~ 4 : 1 and that neu- 
tron stars are even less numerous. A flatter IMF than the 
Miller-Scalo one we have used could increase the ratio of 
white dwarfs to MS stars but we found that such an IMF 
was not well able to reproduce the observed Galactic Cen- 
tre giant population. We do not think it is likely that the 
number of white dwarfs in the Galactic Centre exceeds the 
number of MS stars. It is also possible that a significant 
fraction of neutron stars would be ejected from the Galactic 
Centre by their natal kicks. 

The smallest values of Rmin considered by 

IBailev fc Daviesl (ll999T ) corresponded to ~ -Rgiant/4 for their 
two models. We repeated our three-body/Bondi-Hoyle anal- 
ysis using 1 M impactors and found that collisions in which 
Rmhx < -Rgiant/6 may eject the giants' cores and therefore 
strip the cores of sufficient mass to affect the giants' evo- 
lution in a similar manner to encounters with black holes 
at larger periastrons. Although the critical encounter cross 
section for a MS impactor to have similar effects to a BH 
impactor is then much smaller, the number density of MS 



stars at the Galactic Centre is much larger than that of black 
holes, so that encounters with MS stars may have an effect 
on the giant population comparable to or greater than those 
with black holes. 

We performed SPH simulations to check the results 
of our restricted thre e-body analysis. In agreement with 
IBailev fc Davies] l| 19991 ). we found that encounters at i? m in > 
-Rgiant/4 eject little mass, at most ~ 10% of the enve- 
lope. However, as shown in Figure [T4l the SPH results and 
the three-body results both suggest that single impacts at 
smaller periastrons may be able to expel significant quanti- 
ties of mass, or eject the giant core. We caution, however, 
that we have not performed SPH simulations of encounters 
between giants and MS stars in which the separation be- 
tween the giant core and the centre of mass of the MS was 
less than 4 Mq since the treatment of the MS star as a point 
mass may lead to misleading results at such small separa- 
tions. 

In addition, Figure [3] suggests that giants can expect 
to suffer several tens to of order 100 impacts with MS stars 
at Galactocentric radii < 0.1 pc. We therefore repeated the 
Monte Carlo calculations from the previous section allow- 
ing for multiple impacts with a mixture of MS and BH im- 
pactors. We used the collision probabilities from Figure [3]to 
determine the expected numbers of black hole and main se- 
quence impactors at a given Galactocentric radius. We then 
performed Monte Carlo calculations and calculated the cu- 
mulative effect of all the impactors. In Figure [15] we show 
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Figure 11. Plot of the results of collisions between a lMg gi- 
ant star and a 10 Mg black hole (black squares) at velocities of 
400 (solid line), 800 (short-dashed line) and 1200 (long-dashed 
line) km s . The periastron of the collision is given on the x— 
axis and fraction of the envelope remaining after the collision is 
given on the i/-axis. Black lines with squares are the results from 
SPH calculations. The dotted line is the predicted envelope re- 
maining from the reduced three-body treatment, including the 
sound speed inside the envelope when calculating the capture ra- 
dius, taking = 800 km s" 1 . The horizontal lines represents the 
maximum envelope fraction remaining with which this giant does 
not evolve to become visible in the middle and bright bands at 
the Galactic Centre. 
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Figure 12. Plot of the fraction of 1 Mq giants prevented from 
evolving to become visible in the middle band against the number 
of collisions with 10 Mq black holes, assuming that the mass loss 
as a function of collision periastron is given by the results of SPH 
simulations at velocities of 400 (solid line), 800 (short-dashed 
line) and 1200 (long-dashed line) km s -1 . 
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Figure 13. Fraction of 1 Mq giants prevented from evolving to 
become visible in the middle band as a function of Galactocentric 
radius, assuming that all giants arc the 'standard' model with 
Tcorc = 0.51 (solid line), or assuming giants are randomly dis- 
tributed in age along the giant branch (dashed line). The vertical 
dotted line represents the radius within which the middle band is 
depleted. 
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Figure 14. Plot of the fractional envelope mass remaining after 
collisions between a 1 Mq giant star and a 1 Mq main sequence 
star (black squares) at a velocity of 800 km s -1 . The periastron 
of the collision is given on the a;-axis and fraction of the envelope 
remaining after the collision is given on the y— axis. The dotted line 
is the predicted envelope remaining from the reduced three-body 
treatment, including the sound speed inside the envelope when 
calculating the capture radius. The horizontal lines represents 
the maximum envelope fraction remaining with which this giant 
does not evolve to become visible in the middle and bright bands 
at the Galactic Centre. 
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Figure 15. Fraction of 1 Mq giants prevented from evolving into 
the middle band as a function of Galactocentric radius by col- 
lisions with MS stars (solid line) and BHs (short-dashed line), 
assuming collision rates boosted by a factor of two due to ec- 
centricity of giant's orbit around GC, and using 2— D interpola- 
tion to determine number of impactors, relative velocity and thus 
kill fraction as functions of Galactocentric radius for the BH im- 
pactors and assuming a single velocity of 800km s -1 for the MS 
impactors. The long— dashed line is the fraction of 1 Mq giants 
prevented from evolving into middle band by a mixed population 
of black holes and main sequence stars. 
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Figure 16. Fraction of 2 Mq giants prevented from evolving into 
the bright band as a function of Galactocentric radius by collisions 
with MS stars (solid line) and BHs (short-dashed line), assuming 
collision rates boosted by a factor of two due to eccentricity of gi- 
ant's orbit around GC, and using 2-D interpolation to determine 
number of impactors, relative velocity and thus kill fraction as 
functions of Galactocentric radius for the BH impactors and as- 
suming a single velocity of 800kms — 1 for the MS impactors. The 
long-dashed line is the fraction of 2 Mq giants prevented from 
evolving into the bright band by a mixed population of black 
holes and main sequence stars. 



the results of this analysis for the 1 Mq giant. The figure 
depicts, as a function of Galactocentric radius, the fraction 
of 1 Mq giants prevented from evolving to become visible 
in the middle band, assuming BH impactors acting alone 
(short-dashed line), assuming MS impactors acting alone 
(solid line) and assuming a mixture of BH and MS impacts 
in proportions dictated by the relevant collision probabilities 
(long-dashed line) . Figure [16] depicts the fraction of 2 Mq 
stars prevented from evolving into the bright band. 

We find that the MS stars are approximately as effec- 
tive - somewhat more so in the case of the 1 Mq giant - in 
changing the visible giant population as the BH impactors. 
This result demands some explanation. 

In Figure 1171 we plot for 1 and 2 Mq giants, the ratio 
as a function of Galactocentric radius, of the probability of 
the giant being hit by a main sequence star to the probabil- 
ity of it being hit by a black hole. We see that, in the range 
0.01 pc> r > 0.1 pc, this ratio is « 10 for both giants. In 
Figure [181 we plot the fraction of 1 Mq stars prevented from 
evolving to become visible in the middle band as a function 
of the number of impactors, assuming BH impactors alone 
(dashed line) and MS impactors alone (solid line). We see 
that the factor by which the number of MS impactors must 
exceed the number of BH impactors in order to deplete the 
same fraction of giants is at most ~ 5. It is therefore not sur- 
prising, given the assumptions inherent on our model, that 
MS impactors are as effective in depleting giants as black 
holes. 

The cumulative effect of collisions with black holes and 
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Figure 17. Comparison of the ratio of the probability of being 
struck by a main sequence star to that of being struck by a black 
hole as a function of Galactocentric radius for a 1 Mq giant (solid 
line) and a 2 Mq giant (dashed line). 



Red giant stellar collisions in the Galactic Centre 17 




10 20 30 40 

No. impactors 



Figure 18. Comparison of the fraction of 1 Mq giants prevented 
from evolving into the middle band against the number or im- 
pactors, where each impactor is either a 10 Mq black hole or a 
1 Mq main sequence star. 

main sequence stars on the 1 Mq giants which dominate 
the middle band is to evacuate a three-dimensional volume 
~ 0.04 pc in radius of giants. The effect of collisions on the 
2Mq giants dominating the bright band is smaller, evacu- 
ating a region ~ 0.02 pc in radius. 

7 DISCUSSION 

7.1 Main sequence versus black hole impactors 

We used the results of Monte Carlo modelling combined with 
SPH calculations and three-body simulations to quantify 
the mass lost by a given giant in an encounter at a given 
Rmin with a black hole or a main-sequence star. Taking suc- 
cessive collisions to be independent and assuming that BH 
impactors and MS impactors acted alone or in concert, we 
found that the effects on the visible giant population of the 
two types of impactor were comparable, owing to the much 
greater number density of MS stars than black holes in the 
Galactic Centre. 

This result should be treated with some caution, how- 
ever. In the three-body and SPH calculations, we have 
treated the impactors as point masses. While this is clearly 
valid for the black holes, it may not always be so for the 
MS stars, particularly in the encounters at very small 7? m i n 
required to expel the giant core. 

The periastron necessary to notionally eject the core 
may be so small that it is not valid to approximate the MS 
star as a point mass, either because it is physically too small 
to fit, or because it approaches the core so closely that it ex- 
periences Roche-lobe overflow onto the core. 

The size of the Roche lobe is given by (|EggletonH l983) 

rk = 0.49(m < /m j )§a ^ 

0.6(mi/rrij) 3 + ln(l + (rrn/nij)^) 



Given the core masses given in Table 1, we find that, for 
all the giants, the MS star cannot approach the core closer 
than 2.0 — 2.2 Rq without experiencing Roche lobe overflow. 
The consequences of this phenomenon are not obvious. 
The high-velocity of the intruding MS star may render the 
amount of mass actually transferred in the flyby with the 
giant core very small. In this case, the MS star will exit 
the giant envelope relatively unscathed and the giant core 
will probably be ejected. However, it is possible that the 
tidal distortion of the MS star (or simply an impact at a 
periastron less than the radius of the MS star) will lead to a 
collision between the core and the MS star. This encounter 
may dissipate enough energy that the MS star becomes 
bound to the core and sinks or spirals into the centre of 
the giant to smother the core. The ultimate result of such 
an encounter may be that the MS star becomes part of 
the giant envelope, resulting in a larger and brighter giant. 
Understanding which (if either) of these scenarios is correct 
and what effect they would have on the giant's evolution 
require detailed high-resolution hydrodynamic calculations 
which are unfortunately beyond the scope of this paper, 
but would make an interesting topic of further work. Our 
three-body analysis shows that it becomes significant for 
giants of 3 Mq and above, as the i? m i n required for a 1 Mq 
object to expel core is sa 2Rq. The fractions of giants 
prevented from evolving to become visible in the depleted 
K-bands by MS impactors should therefore be regarded as 
upper limits and it may be that encounters with black holes 
are in reality significantly more important. 



7.2 Encounters with main— sequence binaries 

iDavies et al.l |l998) studied the effects of encounters be- 
tween MS binaries and giants. Since binaries have much 
larger cross sections than single stars, even a relatively small 
population of binaries could be involved in a significant num- 
ber of collisions. The maximum binary cross section that 
can be considered is set by the requirement that the binary 
must be hard (so that it is not disrupted by encounters) 
and is therefore determined by the local velocity dispersion. 
As the velocity dispersion increases, the hard-soft bound- 
ary semimajor axis and therefo re the cross s ection of the 
binary decreases. Additionally, IDavies et al.l (|l998h found 
that, at higher relative velocities, a smaller fraction of inter- 
actions resulted in outcomes d estructive to the gia nt. The 
relative velocities considered by IDavies et all l|l998t) ranged 
from 50 — 150 km s -1 . We extrapolate their results to the 
vel ocities considered in t his p aper of ~ 1000 km s -1 . Follow- 
ing lHut fc Bahcalll l|l985t) and lDavies et all l|l998l) . we define 
the maximum interaction separation s max as a function of 
Voc as 

Smax = ( + 0.6 ) d (17) 

V Voo /VcTit J 

where v cr it is the critical velocity of the chosen binary, 
defined as the velocity at infinity an intruder of the same 
mass as the binary components must have such that the 
total energy of the three-body systems is zero, and d 
is the binary separation. If we ad opt the same binary 
parameters as IDavies et al. | HHl), d = 23.5 R© and 
lie^kms -1 . If Vao is 1000 km s , we obtain the 
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cross section of the binary si rlax ~ 630 R Q . Even if we 
make the assumption that all interactions at i? m i n < s ma x 
prevent the giant evolving into the d epleted K-ba n ds, fo r 
the 93 R0 2 Mq giant considered by iDavies et al.l (| 1998T ) , 
this is « -R|i an t/13. Given that the cross section for a single 
encounter with a 1 Mq impactor to prevent this giant evolv- 
ing into the depleted K-bands is ss -Rgi ant /36, this implies 
that binary interactions are a factor of ~ 3 more effective 
in altering the giant population, implying that the stellar 
population must consist of at least 25% binaries in order 
that binary encounters are more significant th an encounters 
with single MS stars. In agreement with IDavies et al.l 
(1998), we find that the fraction of binaries would have to 
be unrealistically high in the Galactic Centre to have a 
significant collisional effect on the population of giants there. 



7.3 Effect of collisions on the observed giant 
population 

For a single impact with a given giant, more mass has to be 
ejected to prevent the giant evolving into the middle band 
than to prevent it evolving into the bright band, since the 
former entails a more severe perturbation to the star's evo- 
lution. The collision cross section to prevent a giant evolving 
into the bright band is correspondingly larger than that to 
prevent the same giant from evolving into the middle band. 
However, ejecting a given fraction of the envelope from a 
2M0 giant is intrinsically more difficult than ejecting the 
same envelope fraction from a 1 Mq star, since the 2 Mq gi- 
ant's envelope is more strongly bound. In addition, the times 
which stars spend on the giant branch decrease strongly as 
the stellar mass increases. These two facts largely explain 
why we find that the bright band is depleted to a similar ra- 
dius as the middle band. The bright band is dominated by 
2 — 3 Mq giants which have more strongly bound envelopes 
and which are less likely to be struck by MS or BH impactors 
while on the giant branch. 

We have calculated for each model giant the probability 
at a given Galactocentric radius that a combination of col- 
lisions with black holes and MS stars will prevent the giant 
evolving to become visible in the middle or bright bands. 
However, in order to compare our results with observations, 
we must examine the effect of these collisions on the pro- 
jected surface density of sources at the Galactic Centre, since 
this is what is actually observed. If we assume that the gi- 
ant stars are distributed in the same way as the MS stars, 
i.e. that n g i an t(»') oc r~ 1A and that the distribution has a 
maximum radius r max sufficiently large that changing it by 
factors of two does not influence the result significantly, we 
can simply integrate the surface density E of stars at any 
projected radius x as 



E(x) = 



[l-f(z)} n(z)dz, 



(18) 



where z — (r 2 — x 2 ), z n 



x 2 ) and f(z) 



f(r) (where r — (z 2 + x 2 )) is the fraction of giants 
prevented from evolving into the relevant band at three- 
dimensional radius r. The quantity of interest is then the ra- 
tio, as a function of projected radius x, of E (collisions) /E(no 
collisions), since this reveals the effect of collisions on the ob- 
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Figure 19. Surviving fraction of IMq giants visible in the middle 
band plotted against projected Galactocentric radius, assuming 
our standard BH cusp (short— dashed line) and a cusp with four 
times as many BH (long— dashed line). We took r max to be 10 pc. 
Results with r max = 5 pc and r max = 20 pc are very similar. 



served surface density of sources. The results of performing 
this analysis on the 1 Mq giant in the middle band are shown 
in Figure 1191 The two curves show the depletion expected 
due to collisions with MS stars and (i) our standard popu- 
lation of BH (5000 within 0.1 pc, short-dashed line), (ii) a 
black hole population enhanced by a factor of four. We find 
that the standard population of black holes (in combination 
with the MS stars) depletes 50% or more giants from the 
middle band within a projected radius of ~ 0.02 pc, while the 
enhanced BH population depletes 50% or more middle band 
giants within ~ 0.04 pc. This is comparable in size to the de- 
pleted region observed in the middle band bv lGenzel et al.l 
l| 19961 ). The black hole population requir ed to achieve this is 
somew hat larger than that predicted b y [Alexander fc Livid 
l|2004l ) and lDeegan fc Navakshinl (120071 ). 



7.4 Monte Carlo simulations of the Galactic 
Centre 

To examine the influence of collisions on the Galactic 
Centre environment in more detail, we used our Monte 
Carlo population synthesis model of the Galactic Centre 
cluster (from which Figure [2] was constructed). We as- 
sume that star formation there has been proceeding at 
a constant rate for 14Gyr and that stars are born with 
a Miller-Scalo IMF. As discussed in Section 2, when we 
varied these assumptions, we obtained a stellar population 
that was a significantly worse fit to the observed numbers of 
giants in all three K-bands throughout the whole Galactic 
Centre cluster, so we conclude that these assumptions 
are reasonabl e. We use the same STARS evolution tracks 
and Ijohnsor] (| 19661 ) colours and bolometric corrections 
to calculate the evolutionary phase and K-magnitude of 
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Figure 20. Cumulative number of stars against projected Galac- 
tocentric radius measur ed at 0.04, 0,08 and 0.1 2 pc. Squares are 
observed numbers from lTrippe fc Genze j ll2007f) with ■Jn errors. 
Lines are results of Monte Carlo realisations with no collisions 
(solid lines), collisions with our standard black hole and main- 
sequence populations (long-dashed lines) and collisions with a 
black hole population enhanced by a factor of four (dotted lines). 
Thick lines show the the mean numbers of stars averaged over fif- 
teen Monte Carlo realisations and thin lines are +/- one standard 
deviation. 



each star. The stars are distributed in a s phere 2pc in 
radius according to the power laws derived bv lGenzel et al.l 
(2003) using the normalisation described in Equation [3] We 
constructed model clusters neglecting the effects of stellar 
collisions, models including collisions with MS stars and 
our standard and enhanced populations of black holes. We 
then compare d all our models with ne w observational data 
obtained from iTrippe fc Genzell (120071 ). 

We constructed plots of the mean cumulative number 
of stars in our Monte Carlo models in the middle band 
against projected Galactocentric radius and compared them 
to the same plot generated from the observational data. 

In Figure [20] we show the cumulative number of stars 
in the middle band against projected Galactocentic radius 
generated by our Monte Carlo models (lines with +/- la 
standard deviation) and compare these with the observed 
numbers of objects (squares, also with +/- la errors). It 
is clear that both models in which collisions are included 
(dashed and dotted lines in Figure I20[) give a much better 
fit to the observations than a model in which collisions are 
not included (solid line in Figure 120 p , and that the model 
with the enhanced black hole population gives a good fit 
to the observed data. We therefore conclude that collisions 
with MS stars and black holes may be able to account for 
the depletion of giants in the middle band, although with a 
rather large black hole population. 

Repeating this analysis for the bright band (dominated 
by 2 Mq giants) , we find that the region evacuated of 
giants by collisions is too small (~ 0.02 pc in radius) to 
explain the observed depletion, even using our enhanced 
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Figure 21. Plot of total mass densities against Galactocentric 
radius for models with our standard BH population (long dashed 
line), models with black hole populations four (dash-dot line) and 
thirty (short-dashed line) times larger, and after Schodel et al., 
(2007), Figure 7 (solid line). 



population of black holes. We determined that the black 
hole population required to deplete the 2Mq giants out to 
a radius of 0.1 pc is at least thirty times more populous 
than our standard model, so that there would be 1.5 x 10 5 
BH within 0.1 pc. The effect of this population would be to 
deplete the middle band giants out to a radius of ~ 0.3 pc, 
which is much greater than the observed depletion radius. 
This BH population is also much la r ger than that inferred 
by, e. g., lAlexander fc Liviol (|2004j ). Deegan fc Navakshinl 



l|2007h . In addition in Figure 1211 we plot the total mass 
density against Galactocentric radius for the model with 
our standard black hole population and models with black 
hole populations fou r and thirty times la rge, and compare 
it to that inferred in ISchodel et all (|2007h . We see that, at 
0.1 pc, the model with four times our standard number of 
black holes has a m ass density ~ 2 times that inferred by 
ISchodel et all (|2007l ) and the model with thirty times the 
standard number of black holes has a mass density ~ 10 
times greater. We therefore conclude that the observed 
depletion of the brightest giants in the Galactic Centre 
cannot be due purely to collisions of giants with MS stars 
or black holes. The depletion of the brightest giants may 
instead be due to tidal stripping of stars on very eccentric 
orbits, and it is also possible that modification of the stellar 
population by collisions of MS stars with each other may 
decrease the numbers of low- and intermediate-mass giants 
that contribute to the bright K-band. We will look at these 
scenarios in detail in subsequent papers. 

We also checked to see what effect collisions would 
have on giants in the undepleted faint band (15 >K> 12). 
This band is dominated by 1 — 2M© objects. Preventing 
giants of these masses evolving into this brightness band is 
extremely difficult, requiring the loss of > 99% of the giant's 
envelopes - even with the enhanced black hole population 
of 2 x 10 4 within 0.1 pc, the radius to which faint band 



20 James E. Dale, Melvyn B. Davies, Ross P. Church, Marc Freitag 



objects would be depleted is < 0.01 pc. We conclude that 
collisions with black holes and MS stars are unlikely to 
have an observable effect on the fainter population of giants. 



8 CONCLUSIONS 

Our main conclusions may be summarised as follows: 

(i) The stars missing from observations in the bands 
K < 10.5 and 10.5 < K < 12 in the Galactic Centre are 
low-mass giant stars, in the mass range 1 — 4Mq. 

(ii) The Galactic Centre is a highly collisional environ- 
ment. In particular, within a Galactocentric radius of 0.1 pc, 
giant stars are likely to suffer multiple collisions with both 
main-sequence stars and black holes. 

(iii) We have shown that significantly altering the evolu- 
tion of giant stars by means of collisions requires the ejection 
of large quantities of mass (> 20% of the envelope) while the 
giant is on the first giant branch. Removal of ~ 20% of the 
envelope of a IMq giant is sufficient to make the AGB sig- 
nificantly fainter, but has little effect on the RGB phase. 
Mass losses of > 40% occurring on the RGB are sufficient to 
decrease the brightness of the RGB tip by ~ 1 magnitude 
and prevent the star evolving onto the HB or AGB. Mass 
loss on the horizontal branch and AGB has essentially no 
effect on a star's subsequent evolution, as we showed in Fig- 
ure [5] 

(iv) We have shown that penetrating encounters with 
10 Mq black holes expel the cores of giants. The cores carry 
away some fraction of the envelope, resulting in very large 
fractional mass losses. Such mass losses in single collisions 
can affect the evolution of the giants to such a degree that 
they never evolve to become brighter than K = 12. 

(v) We have also found, using a restricted three-body 
analysis, that encounters at very small i? m i n with 1 Mq 
main-sequence impactors can also eject the giant core, lead- 
ing to the severe mass losses required to prevent the giant 
evolving to become visible in the depleted K-bands. 

(vi) The cumulative effects of collisions with black hole 
and main-sequence impactors are able (for a population of 
black holes of 2 x 10 4 within 0.1 pc) to deplete a region 
« 0.04 pc (« 1 arcsec) in radius of giants in the middle K- 
band, comparable in size to the region observed to be devoid 
of these objects bv lGenzel et al.l (|l996l ). Collisions can plau- 
sibly explain the depletion in this band. 

(vii) The cumulative effects of collisions deplete the 
2 — 3 M© giants which dominate the bright K-band to a 
somewhat smaller radius (~ 0.02 pc). This is considerably 
less than the si ze of the region obs erved to be depleted of 
these giants by iGenzel et al.l (| 1996h . The depletion of the 
very brightest giants cannot be explained by collisions with- 
out invoking a black hole population so large that it would 
deplete the middle band out to radius much larger than ob- 
served. 

(viii) Owing to the extreme mass losses required to pre- 
vent giants evolving into this band, collisions are not likely 
to have an observable effect on the faint band (15 >K> 12). 
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